Breathers and surface modes in oscillator chains with Hertzian interactions 
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Abstract — We study localized waves in chains of oscil- 
lators coupled by Hertzian interactions and trapped in local 
potentials. This problem is originally motivated by New- 
ton's cradle, a mechanical system consisting of a chain of 
touching beads subject to gravity and attached to inelastic 
strings. We consider an unusual setting with local oscil- 
lations and collisions acting on similar time scales, a sit- 
uation corresponding e.g. to a modified Newton's cradle 
with beads mounted on stiff' cantilevers. Such systems sup- 
port static and traveling breathers with unusual properties, 
including double exponential spatial decay, almost van- 
ishing Peierls-Nabarro barrier and spontaneous direction - 
reversing motion. We prove analytically the existence of 
surface modes and static breathers for anharmonic on-site 
potentials and weak Hertzian interactions. 

Granular media are known to display a rich dynami- 
cal behavior originating from their complex spatial struc- 
ture and diff'erent sources of nonlinearity (Hertzian con- 
tact interactions between grains, friction, plasticity). In 
the case of granular crystals (i.e. for grains organized on 
a lattice), nonlinear contact interactions lead to diff'erent 
types of localized wave phenomena that could be poten- 
tially used for the design of smart materials such as acous- 
tic diodes Q. Among the most studied types of excita- 
tions, solitary waves can be easily generated by an impact 
at one end of a chain of touching beads (see |7 , 8 1 and ref- 
erences therein). In the absence of an original compression 
in the chain (the so-called precompression), these solitary 
waves diff'er from the classical KdV-type solitary waves, 
since they are highly-localized (with super-exponential de- 
cay) and their width remains unchanged with amplitude 
(see e.g. 1 12|). These properties originate from the fully 
nonlinear character of the Hertzian interaction potential 
V(r) = f r(-r)+/^ (with y > and (a)+ = max(a,0)), 
which yields a vanishing sound velocity in the absence of 
precompression. 

Discrete breathers (i.e. intrinsic localized modes) form 
another interesting class of excitations consisting of time- 
periodic and spatially localized waveforms [6, 3|. These 



waves exist in diatomic granular chains under precom- 
pression ||2| [131 (with their frequency lying between the 
acoustic and optic phonon bands) and can be generated 
e.g. through modulational instabilities. However, because 
precompression suppresses the fully nonlinear character of 
Hertzian interactions, these excitations inherit the usual 
properties of discrete breathers, i.e. their spatial decay is 
exponential and their width diverges at vanishing amplitude 
(for frequencies close to the bottom of the optic band). 

The situation is sensibly diff'erent for granular systems 
without precompression. In that case, localized oscillations 
can be generated on short transients in the form of transi- 
tory defect modes induced by a mass impurity (see [ llj 
and references therein) but never occur time-periodically 
as proved in |5 1. Indeed, uncompressed granular chains are 
described by the Fermi-Pasta-Ulam lattice with Hertzian 
interactions 



% = V'(yn+i - yn) - V'iyn - yn-i) 



(1) 



(or spatially inhomogeneous variants thereof), where yn(t) 
denotes the ^th bead displacement from its reference posi- 
tion. For all T-periodic solutions of ([B, the average inter- 
action forces V'(yn+i -Jrddt are independent of n (this 
is immediate by integrating ([B). Consequently, localized 
oscillations would yield a vanishing average interaction 
force between grains, which is impossible since Hertzian 
interactions are repulsive under contact and vanish other- 
wise. 

In contrast to the above picture, we have numerically 
established in | 5 1 the existence of time-periodic localized 
oscillations in Hertzian chains with symmetric local poten- 
tials described by the system 

% + W\yn) = V\yn^l - yn) - y\yn - yn-l) (2) 

where W{y) = | + | and 5- g R measures the local an- 
harmonicity. System (|2]) describes small amplitude waves 
in a Newton's cradle |4 1 (figure[T] left) or other mechanical 
systems consisting of beads mounted e.g. on an elastic ma- 
trix [10] or cantilevers [5] (figure [T] right). In the last two 



cases, local oscillations and collisions between beads can 
occur on similar time scales for realistic material parameter 
values, which allows for breathing dynamics to take place. 
Static and moving breathers can be generated from stan- 
dard initial conditions such as a localized impact 1 5 , JOJ or 
perturbations of unstable periodic traveling waves (4]. 
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Figure 1: Left : prototypical Newton's cradle consisting 
of a chain of beads attached to pendula. Right : array of 
clamped cantilevers decorated by spherical beads. 

Let us summarize the results of [5|. Two fami- 
lies of static breather solutions of ^ (parametrized by 
their frequency ojb > I) have been computed using the 
Gauss-Newton method and path-following. The first one 
consists of bond-centered breathers satisfying y^it) = 
-y-n+i(t) := Siyn(t), and the second corresponds to site- 
centered breathers with y^it) = -y-n(t + Tiy/2) := S2yn0), 
where denotes the breather period (each solution fam- 
ily is invariant by a symmetry Si of Profiles of both 
breather types are given in figure |2] for a particular fre- 
quency close to unity (i.e. the linear frequency of lo- 
cal oscillators), a limit where the breather amplitude van- 
ishes. Near this limit, quasi-continuum approximations of 
the breather profiles can be derived [5 1, namely 

/^\t) = 2€(-ir[g(n) + g(n-l)] cosM (3) 

for bond-centered breathers and 



breather), for anharmonic on-site potentials, but more criti- 
cally even in the harmonic case. This phenomenon is linked 
with an extremely small diff'erence between the energies 
of site- and bond-centered breathers having the same 
frequency (the so-called approximate Peierls-Nabarro bar- 
rier). Typical values of breather energies are given in figure 
[3 where *H = 2« ^« and en = ^yl-\- W(yn) + V(yn+i - yj- 
Another manifestation of the high breather mobility is the 
systematic formation of a travelling breather after an im- 
pact at one end of the oscillator chain (see figure |4]). 
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Figure 2: Bond-centered (left plot) and site-centered (right 
plot) breather solutions of ^ for oj^ = 1.01 (potential pa- 
rameters are 7 = 1 and ^ = 0). The numerical solutions 
(marks) are compared to the quasi-continuum approxima- 
tions yn\ y^n^ (continuous lines). 




yf(f) = 2e i-ir [^{n + ^) + g(n - ^)] cos (Wb?) (4) 



for site-centered ones, where aj\, 

x3 



g(x) = (^y cos^(f)for|x| < f and^ 



1 + 1;^ (TO ^ 1.545), 
elsewhere. 



Figure |2] compares the numerical solutions with the above 
approximations for oj^ ^ I, showing excellent agreement. 
Obviously discrepancies appear at larger amplitude, since 
approximations (O and retain only one Fourier mode 
and are independent of the local anharmonicity (which is 
dominated by the Hertzian nonlinearity at small ampli- 
tude). The above approximations possess compact supports 
(with a width independent of ojh), but the breathers com- 
puted numerically do not share this property strictly speak- 
ing and instead display super-exponential localization (this 
is in analogy with the case of homogeneous polynomial in- 
teraction potentials, see section 4.1.3 of |3 |). 

Dynamical simulations of iQ indicate that extremely 
small perturbations of the static breathers can lead to 
their translational motion (generating a so-called traveling 



Figure 3: Breather energies versus frequency ojt, when the 
on-site potential is harmonic (s = 0, left plot) or anhar- 
monic (s = -1/6, right plot) and 7=1. The red curves 
give the energy "Hbc of bond-centered breathers, the blue 
curves the approximate Peierls-Nabarro barrier Epn (en- 
ergy diff'erence bewteen site- and bond-centered breathers) 
and the black curves the relative energy ratio EpN^Htc- For 
small amplitude breathers (i.e. ojh ^ I), the diff'erent values 
of s yield comparable values of Ep^. Clearly Epn increases 
with breather amplitude but remains very small in this pa- 
rameter range (e.g. Ep^ is close to 10"^ for coh = 1.5 and 
s = -1/6). The harmonic case yields even smaller barriers, 
by 3 - 4 orders of magnitude for cl)^ = 1.3. 

Diff'erent phenomena have been identified in |5 1 depend- 
ing on the softening or hardening character of the local 
potential W. Firstly, the stability of both site-centered 
and bond-centered breathers is critically dependent on the 
strength (and sign) of the anharmonicity (we refer to [5 1 for 
more details). Secondly, depending whether ^ < or ^ > 0, 



the occurence of surface mode excitations (i.e. oscillations 
localized near a boundary) or direction-reversing traveling 
breathers was observed after an impact. Both phenomena 
are illustrated by figures |4] and |5] The origin of direction 
reversal is still unclear at the present stage, although we 
think it might originate from the interaction between the 
traveling breather and nonlinear waves confined between 
the breather and the boundary (figure |5] right plot). 




Figure 4: Space-time diagrams showing the interaction 
forces fn = -(jn - yn+\f'+^ in system Q with 7=1 and 
free end boundary conditions. Forces are represented in 
grey levels, white corresponding to vanishing interactions 
(beads not in contact) and black to a minimal negative value 
of the contact force. The initial condition is ^^^(O) = for 
n > 1, jniQ) = for /2 > 2, ji(0) = 0.94. For an harmonic 
local potential {s = 0, left plot), a traveling breather is gen- 
erated. For a soft local potential (s = -0.7, right plot), one 
observes in addition the excitation of a surface mode. 




Figure 5: Left : space-time diagram of the energy density 
en for system ^ with y = I and a hard anharmonic local 
potential (s = 1). The initial condition is the same as in 
figure m except ji(0) = 1.9. Right : corresponding particle 
velocities close to direction-reversing (t 241). 

In order to understand the absence of surface mode ex- 
citation observed in figure |5] for ^ = 1 , we now analyze 
the properties of surface modes computed by the Gauss- 
Newton method (all numerical computations are performed 
for 7 = 1). The results are shown in figure |6] (left column). 
For ^ = 1, surface modes exist for frequencies ojs lying 
above a critical frequency ajj^m ~ 1-96, where a pair of 
Floquet eigenvalues converges towards unity (spectra not 
shown here). For ojs ~ o^min these solutions are spectrally 
stable and one can observe an energy threshold. This ex- 
plains why a surface mode is not observed in figure|5]where 



the energy of the initial excitation is well below the excita- 
tion threshold. This situation contrasts with the case of soft 
local potentials illustrated in figure |6] (right column). For 
s = -0.7, surface modes exist for ojs e [0.705, 1). They 
are spectrally stable for ojs ^ I but oscillatory instabili- 
ties appear at smaller frequencies (spectra not shown here) . 
When ojs ^ I, their energy and amplitude vanish, opening 
the possibility of exciting such modes for arbitrarily small 
initial velocities of the first particle. 




Figure 6: Top : profile of a surface mode at the instant 
of maximal amplitude, for a hard local potential (s = 1, 
ojs = 2, left plot) and a soft one (s = -0.7, ojs = 0.99, 
right plot). Bottom : energy of the surface mode versus 
frequency oJs, for ^ = 1 (left plot) and s = -0.7 (right plot). 

Having numerically analyzed the properties of surface 
modes of (|2]) when y (the stiff'ness constant of Hertzian 
interactions) equals unity, we now prove the existence of 
time-periodic and spatially localized solutions when y is 
small. The small coupling (anticontinuum) limit was in- 
troduced in reference [6 | (and considerably generalized in 
1 9 1) and applies both to breathers and surface modes. From 
a physical point of view, this parameter regime can be re- 
alized with an array of clamped cantilevers (see figure [T) 
decorated by spherical beads made from a sufficiently soft 
material (e.g. rubber). In addition, the anticontinuum limit 
requires an anharmonic local potential W, i.e. s ^ 0. 

We consider system Q for ^ g F, with F = Z for discrete 
breathers (case of a doubly infinite chain) and F = Mo for 
surface modes (semi-infinite chain, with free-end bound- 
ary conditions at ^ = 0). For 7 = 0, system Q admits 
localized periodic solutions satisfying = ^ 0, 
yo + W\yo) = and yo(t + Tq) = yo(t). Setting yo(0) = 0, 
a := yo(0) > and Tq = To(a), it is a classical result that 
s T^ia) < 0. Consequently, the time periodic oscillations 
yo can be parametrized by their period Tq, with Tq < In 
for ^ > and To > 2;: for ^ < 0. We shall denote by 
y^jQ '•= (yn)n the even To-periodic solution constructed 



above for y = 0. Now let us assume Tq ^ InM (this 
condition is automatically satisfied for s > 0). Thanks to 
this nonresonance condition, one can apply the persistence 
theorem for normally non-degenerate discrete breathers 
proved in reference [9 1 (section 3.5), which readily applies 
to system Q since V e C^(R) has sufficient smoothness. 
More precisely, for y small enough in system (|2]), this re- 
sult ensures the existence of a family of localized solutions 
Yyj(. -h 0) parametrized by their period T ^ Tq (since 
9^ 0) and an arbitrary phase 0. In addition the per- 
turbed and unperturbed solutions are close when y is small, 
i.e. one has Yyj = Yqj -h 0(y) in C^(S^'j2(r)) (we de- 
note S ^ = R/(TZ)). Note that the corresponding breather 
solutions are site-centered, due to the local uniqueness of 
y Yyj and the invariance of (|2]) by the symmetry 82- 
Similarly, starting from a bond-centered breather for 7 = 
(with yi = -yo, y^ = elsewhere) allows one to obtain 
bond-centered breathers for 7 ^ 0. More generally, multi- 
breathers having many excited sites can be constructed in 
the same way at small coupling. 

Now let us estimate the spatial decay of breathers and 
surface modes. Consider a T-periodic solution of Q sat- 
isfying the nonresonance condition T ^ 2nM. We assume 
this solution is spatially localized, i.e. its supremum norm 
Vn = WynWoo := sup^^[0 7^j |j„(OI satisfics lim|;,|^+oo v„ = 0. In 
order to estimate v„, we rewrite (|2]) in the form 

yn=L-' [y(y^_,-y^)l^^-y(y^-y^^,)l^^-Syl), (5) 

where L = dj I : C^(S^;R) C^(S^;R) is invertible 
thanks to the nonresonance assumption. By using (|5]), there 
exists C > and m g Mo such that for |^| > m -h 1 we 
have Vn < ^ + ^n^^ + ^H-i)- Introducing the sequence 
€j = sup|^|>y Vyt, it follows that €j < Cej^^ for 7 > m -h 1, 
hence one obtains by induction 

6,<C-2i(3/2y, i = (c2.„p", j>m. (6) 

Since lim^^+oo ej - 0, we can fix m such that X <\. Esti- 
mate © implies ||j«||co < e\n\ < ^(3/2)i"i |^| ^ 
the solution decays doubly exponentially at infinity. 

As a conclusion, we have reviewed recent results of 
||5l on the numerical observation of breathers and sur- 
face modes in granular chains with stiff local potentials. 
We have found the existence of an energy threshold for 
hard surface modes, which explains their non-excitation by 
moderate impacts reported in [5 1. We have obtained dou- 
bly exponential decay estimates for breathers and surface 
modes and proved their existence for anharmonic local po- 
tentials when y ^ 0. The persistence of these solutions for 
larger values of y is numerically observed (for 7=1) but 
remains an open question from an analytical point of view. 
In addition, the spontaneous direction-reversing motion of 
discrete breathers occuring for hard local potentials and the 
smallness of their Peierls-Nabarro energy barrier has yet to 
be understood. 
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